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Abstract 

Recent work on the dynamics of a crystal surface [T. Frisch and A. Verga, Phys. 
Rev. Lett. 96, 166104 (2006)] has focused the attention on the conserved Kuramoto- 
Sivashinsky (CKS) equation: dtu = —dxxiu + Uxx + Ux), which displays coarsening. 
For a quantitative and qualitative understanding of the dynamics, the analysis of 
steady states is particularly relevant. In this paper we provide a detailed study of the 
stationary solutions and their explicit form is given. Periodic configurations form an 
increasing branch in the space wavelength-amplitude (A-^), with dX/dA > 0. For 
large wavelength, A = A^/A and the orbits in phase space tend to a separatrix, which 
is a parabola. Steady states are found up to an additive constant a, which is set by 
the dynamics through the conservation law dt{u{x,t)) = 0: a(A(t)) = A^(t)/48. 

Key words: Nonlinear dynamics. Coarsening, Instabilities, Crystal growth 
PACS: 05.45.-a, 45.20.D-, 68.35.Ct, 81.10.Aj 



1 Introduction 



The study of the growth processes of crystal surfaces [1,2,3] has turned out to 
be source of a variety of nonlinear dynamics. A first, general distinction should 
be made between a crystal growing along a high symmetry orientation (e.g., 
the face (100) of silicon) and one growing along a vicinal orientation (e.g., the 
face (119) of copper). In the former case, growth proceeds [4] via nucleation, 
aggregation of diffusing adatoms and coalescence of islands. The growth is 
intrinsically two-dimensional and the variable describing the growth dynamics 
is the local height z{x,y,t) of the surface, or alternatively, the local slope 
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m(x, y, t) = Vz{x, y, t). The ideal growth of a high-symmetry surface is called 
layer-by-layer. In the latter case (vicinal orientation), the surface is made up 
of a train of steps [5] which capture diffusing adatoms and advance. Since 
steps are neither created nor destroyed, a vicinal surface is usually described 
in terms of step positions Q{x,t), where i labels the steps. In this case the 
ideal growth is called step-flow. 

The interest in the nonlinear dynamics of a crystal surface comes from the 
observation that growth is often unstable [6]. Layer-by- layer growth may be 
destabilized by mound formation, whose phenomenology has many similari- 
ties with thermal faceting [7], even if the driving force is kinetic rather than 
thermodynamic. The representative equation, 

dtZ^-V-[V^m + j(m)] , (1) 

has strict resemblances with the Cahn-HiUiard equation [8] , once we take the 
gradient of both sides. 

The growth of a vicinal surface allows for more rigorous treatments and the 
original two-dimensional character of the growth can be often decomposed in 
the two spatial degrees of freedom, the coordinate x of the single step and the 
index i for different steps, according to the type of emerging instability. In 
fact, step flow may undergo two different types of instability: step bunching 
and step meandering. Step bunching means that the density of steps docs not 
keep constant because steps prefer to gather in bunches separated by large 
terraces. Step meandering means that steps do not stay straight and start 
wandering. 

The possibility that both instabilities appear at the same time is very rare [9] : 
typically, step bunching occurs while steps are straight and step meandering 
occurs while different steps wander in phase and their distance remains con- 
stant [10]. This means that the dynamics can be described by effective one- 
dimensional equations. For example, step bunching in the presence of large 
desorption can be described by the Benney equation [11], 

dfU — tJ'XX ~l~ ^"^xxx fJ'xxxx ~l~ '^x 1 (^) 

giving chaotical dynamics or a regular structure of stable bunches, depending 
on the strength of the propagative 7— term. 

The dynamics of meandering depends on the asymmetry in the attachment 
kinetics to the steps: meandering occurs when adatoms preferentially attach 
to the ascending steps [12]. In the presence of strong evaporation, we get the 
well known Kuramoto-Sivashinsky equation [13], corresponding to the Benney 
equation when 7 = and producing spatio-temporal chaos. In the opposite 
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limit of vanishing evaporation and large asymmetry, a strongly nonlinear gen- 
eralization of the Cahn-Hilliard equation was found [14] : 

dfU = -A{u)dxx[B{u) + C{u)uxx\ ■ (3) 



This equation may produce [15] a constant-wavelength pattern (A = Aq) with 
diverging amplitude or a coarsening process, i.e., an increase in time of the 
the wavelength A. This increase may be either perpetual, therefore defining 
the coarsening exponent n, X{t) ~ i", or interrupted [16] at some length X^^. 

The case of vanishing desorption and weak asymmetry has been recently con- 
sidered by T. Frisch and A. Verga [17], who found the equation 

dtu = -d^xiu + u^x + ul) , (4) 



which has been called the conserved Kuramoto-Sivashinsky (CKS) equation. 
The equation has been solved numerically and appears to give rise to a coarsen- 
ing pattern, whose exponent, according to similarity arguments [17], is n = |. 
A similar equation, with an extra propagative term jUxxx, arises in step bunch- 
ing dynamics with vanishing desorption [18] and in a completely different do- 
main, sand-ripple dynamics [19]. The propagative term breaks the x —x 
symmetry, but does not seem to change the long-time behavior: perpetual 
coarsening with n = | is preserved. Additional discussion about this equation 
can be found in the last Section. 

Most times, coarsening occurs because the Partial Differential Equation (PDE) 
describing the dynamics has a branch of steady states whose wavelength is an 
increasing function of their amplitude. These steady states are unstable with 
respect to phase fluctuations [20] and the profile evolves in time keeping close 
to the stationary branch. 

In Ref. [17] the authors remark that the profile emerging from the CKS equa- 
tion can be thought of as a superposition of parabolas and that a stationary 

parabola is a particular solution of the equation. In this paper we aim at a 
detailed characterization of the steady states of the CKS equation. This will 
be mainly done analytically, giving the general exact relation between u and 
Ux, i.e., determining the trajectories in phase space. 



3 



2 Steady states 



The steady states of the CKS equation (4) , satisfy the second-order nonhnear 
differential equation 

u + Uxx + = a + 6a; . (5) 



Since the constant h must vanish in order to get hmited sohitions, while the 
constant a can be trivially absorbed into a uniform shift of u{x) and be set to 
zero, the problem reduces to solving the differential equation 

Uxx ^-u-ul . (6) 



Interpreting time, this equation describes the dynamics of a harmonic 

oscillator subject to an additional force proportional to the square velocity It 
can be easily solved by evaluating the derivative of u^x = 



Fx = -Ux - '2UxUxx = -Ux{l + 2F) , 



(7) 



so that 

du 
dF 



1 + 2F ' 



(8) 



whose integration starting from a zero- velocity initial condition, u{0) — A and 
Ux{0) — 0, that implies F{0) — —A, gives the first integral 



2{A -u)=ln 



1 + 2F 



1-2A 



(9) 



or 



1 + 2F = ±(1 - 2A) e^(^-") ; 



(10) 



a continuity argument allows us to drop the minus sign 
F^-u-ul^-^ + {l- 2A) e2(^-") , 



fill 



making it apparent that ii A > ^ the force is strictly negative and the tra- 
jectory is not limited. On the other hand, if A = A^, with < A^ < | , 
the particle initially acquires a negative velocity, but the increasing positive 
contribution to the force F will eventually restore Ux — at some position 
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Fig. 1. Trajectories in the phase space (n, n,^.) of the limiting parabola = ^ — u 
(thick full line) and of the periodic steady states: A_ = 1 (dotted), A_ = 2 (short 
dashed), yl_ = 4 (dot dashed), and yl_ = 8 (long dashed). The thin full lines 
correspond to the approximation = 2A-{u + A_), valid close to u = —A^. 

u — —A^ < 0; the converse occurs starting from A — —A^ < 0, so that a 
periodic orbit occurs for any A < \. The pairs of turning points —A^) 
are solution of 

(1 - e-(i-2^+) = (1 + 2A_) e-(i+2^-) , (12) 



that maps the interval [0, |) onto the interval (— oo,0]. Eventually, we can 
write the trajectories in the phase space {u, Ux) as 

ul^\-u-{\-A) e^('^-) ; (13) 

these are plotted in Fig. 1 for a choice of values of A_. If A^ = |, we obtain 
the separatrix u^. = ^ — u, which corresponds to the parabolic trajectory 
u{x) = I - (x - xoY/4:. 

Eq. (13) can be rewritten in a more familiar form, 

^ + V{u)^E^l, (14) 
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Fig. 2. Effective potentials V{u) for bound states, A < ^ (dotted lines), for the 
separatrix, A = ^ (dash-dotted line), and for unbounded states, A > ^ (dashed 
line). The "energy" of the particle is E = ^ (full line) in all cases. The condition 
V{u) = E identifies the turning points. Inset: the potential V*{h) = + e~'^^) 
(full line) and the energy = \{\- A)-\\n{\- A) (dotted line), valid for ^ < \. 

showing that the system satisfies the conservation of a pseudo-energy E = 
with a potential 

V{u) = '-^ + '^{l-A)e'^^--^ (15) 

that depends on the boundary (initial) conditions. The above discussion about 
closed and open orbits becomes a straightforward analysis of the potential 
shape, as reported in Fig. 2. 

It appears from the figure that the potential keeps its shape: this can be 
checked analytically by considering a shift in the u variable, u{x) — h{x) + qo', 
ioT qo — A + ^ln\-^ — A\ the first integral (14) can be rewritten as 

l!^ + V*{h)^ E*{A), (16) 

with E*{^A) = Id — A) — |ln|| — A\, and the A-independent potential 

V*{h) = ^{h±e-^'') , (17) 
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Fig. 3. Periodic steady states configurations, corresponding to = 1 (dotted), 
A- = 2 (short dashed) , A- = 4 (dot dashed) , and A_ = 8 (long dashed) . The thick 
full lines corresponds to the limiting parabola u = ^ — Small inset: the minima 
of the configurations and the approximations u = —A- + A-x'^/2 (thin full lines). 
Large inset: wavelength A as a function of the amplitude The full line is the 
asymptotic relation A = 4:^/A^. 

where the sign + (— ) is valid ior A < ^ (A > For bounded states, the 
energy mininnim is E*{0) = |(1 + ln2) and coincides with the minimum of 
the potential V*{h=ln2/2) (trivial orbit). The potential V*{h) and a repre- 
sentative value of the energy E* are plotted in the inset of Fig. 2. 

The numerical results for the periodic configurations u{x) in real space are 
shown in Fig. 3. It is useful to determine the amplitude A_ in the negative 
X direction, as a function of A+, through the condition (12). For A+ 0, it 
is easily found that A^/A^ — > 1, while in the important limit A^ — > |, A- 
diverges logarithmically according to 

A_e-'^- ^ , (18) 

i.e., 

A_ -|ln(l-2A+) . (19) 

In proximity of m = —A_, bounded trajectories have a minimum with a cur- 
vature that diverges as A^ ^ oo, as shown by the expansion u — — A_ -|- £ in 
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Eq. (13), which gives 

u ^ -A_ + —{6xf + 0[A^_{6x)^] . (20) 
2 



This approximation is shown as thin full lines, both in Fig. 1 and in Fig. 3. 
The quadratic and quartic terms are of the same order when |5a:| ^ l/yC^T, 
which sets the size of the high-curvature region. In fact, the slope Ux{Sx = 
1/y/AZ) PS vC5r joins to the same slope of the limiting parabola u — | — a;^/4, 
when u ~ —A^. 

As for the unbounded states, A > i, it is worth noting that they diverge 
at a finite value Xoo] in the mechanical analogy, this means that the particle 
escapes to infinity in a finite time. We can determine integrating Eq. (13) 
once more: 

= _ 'f* = 1 / , (21) 



A 



and, for large A, 

1 r ds 71 



(22) 







Finally, in the large inset of Fig. 3 we plot the wavelength A of the steady 
states as a function of their amplitude A_. The full line, A = 4:^/A_, gives the 
analytical approximation valid for large A^. It can be determined from the 
asymptotic parabola m = | — x^/A imposing u{\/2) = — or equivalently, 
from the mechanical analogy for the potential V* (see inset of Fig. 2). 



3 Steady states and dynamics 



In the Introduction we have argued that steady states are important because 
dynamics proceeds evolving along the family of steady states of increasing 
wavelength A. In the case of the CKS equation special attention should be 
paid to the constant a appearing in Eq. (5) and to the conserved character of 
Eq. (4). We are now going to show that the conservation law fixes the value 
of a, as a function of A. Afterwards, we explain that for the dynamics a is the 
time dependent vertical shifting of the surface profile. 

On the one hand, the conserved dynamics requires that the spatial average 
{u{x,t)) is time independent; on the other hand, {u{x)) ^ and does depend 
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on A. Therefore, the family of steady states which is relevant for the dynamics 
is 

uxix) = u{x) + a(A) , (23) 



where the constant a satisfies the condition 

a(A) = -{u{x)). (24) 

In the above notation, u{x) means the general steady state (of wavelength 
A) found in the preceding Section. For large A, the average value of u{x) can 
be safely determined by approximating it with the arc of the (separatrix) 
parabola, so that 

A/2 2 2 

{u{x)) = j J dx(^-^)+ o{X') = + o(A^) . (25) 

-A/2 



Therefore, we get 

«(A) = ^ . (26) 



4 Conclusions and discussion 



The starting and only assumption of this work is that steady states of increas- 
ing wavelength play a central role in the dynamics of the CKS equation and, 
in particular, in the coarsening process. This assumption is based on previous 
work [20] on a plethora of different PDE for which this connection has been 
well estabhshed. It is also supported by numerical work [17] on the CKS eq., 
showing that the system indeed coarsens with a profile made up of a sequence 
of arcs of parabola. 

In Section 2 we have analytically found all the steady states and in Section 3 
we have determined the relation between vertical shifting of the profile and 
coarsening process. The theory of the phase diffusion equation [20], which 
would allow a more rigorous derivation of the coarsening exponent n = | [17], 
can not easily be applied to the CKS equation, because the linear operator 
coming from the Prechet derivative is not self-adjoint, a feature that the CKS 
equation shares with the Kuramoto-Sivashinsky equation. 

The fact that the surface profile is made up of arcs of parabola joined by van- 
ishing regions of diverging amplitude (angular points) suggests an alternative 
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dynamical description: the system is made up of a sequence of "particles" (the 
angular points) which tend to annihilate, and therefore the system tends to 
coarsen. Is it possible to study the effective dynamics of these particles? 

Let us finally comment on a modified version of the CKS equation [18,19], 

dtu = -d^^{u + u^^ + ul + ^u^) , (27) 

where a third-order dispersive drift term has been added to the right-hand- 
side. This equation has been mentioned in the Introduction as emerging in 
two different domains: growth of a vicinal surface subject to a step-bunching 
instability and sand ripple dynamics. According to numerics [18], Eq. (27) 
displays a surface profile which is not much different from the profiles shown 
by the CKS equation. Also, the system seems to coarsen with the same law, 
A ~ 

It may be tempting to treat Eq. (27) along the same lines followed for the 
CKS equation. The first remark is that periodic steady states should now be 
replaced by travelling periodic configurations. In fact, Eq. (27) does not have 
periodic stationary solutions because of the 7— term. Nonetheless, one might 
look for periodic solutions depending on {x — vt). It would be interesting to 
understand if there is a branch of travelling periodic configurations in the 
space X — A and if dynamics coarsens following this branch, as for the CKS 
equation. 
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